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ABSTRACT 

The Geosynchronous Imaging Fourier Transform Spectrometer (GIFTS) Sensor Module (SM) Engineering Demon- 
stration Unit (EDU) is a high resolution spectral imager designed to measure infrared (IR) radiances using a 
Fourier transform spectrometer (FTS). The GIFTS instrument employs three focal plane arrays (FPAs), which 
gather measurements across the long- wave IR (LWIR), short/mid-wave IR (SMWIR), and visible spectral bands. 
The raw interferogram measurements are radiometrically and spectrally calibrated to produce radiance spectra, 
which are further processed to obtain atmospheric profiles via retrieval algorithms. This paper describes the 
GIFTS SM EDU Level IB algorithms involved in the calibration. The GIFTS Level IB calibration procedures 
can be subdivided into four blocks. In the first block, the measured raw interferograms are first corrected for 
the detector nonlinearity distortion, followed by the complex filtering and decimation procedure. In the second 
block, a phase correction algorithm is applied to the filtered and decimated complex interferograms. The result- 
ing imaginary part of the spectrum contains only the noise component of the uncorrected spectrum. Additional 
random noise reduction can be accomplished by applying a spectral smoothing routine to the phase-corrected 
spectrum. The phase correction and spectral smoothing operations are performed on a set of interferogram 
scans for both ambient and hot blackbody references. To continue with the calibration, we compute the spectral 
responsivity based on the previous results, from which, the calibrated ambient blackbody (ABB), hot blackbody 
(HBB), and scene spectra can be obtained. We now can estimate the noise equivalent spectral radiance (NESR) 
from the calibrated ABB and HBB spectra. The correction schemes that compensate for the fore-optics offsets 
and off-axis effects are also implemented. In the third block, we developed an efficient method of generating 
pixel performance assessments. In addition, a random pixel selection scheme is designed based on the pixel 
performance evaluation. Finally, in the fourth block, the single pixel algorithms are applied to the entire FPA. 

Keywords: Remote sensing; Fourier transform spectrometry; Geosynchronous Imaging Fourier Transform Spec- 
trometer (GIFTS); interferogram processing; radiometric calibration; Level IB algorithms. 

1. INTRODUCTION 

The Geosynchronous Imaging Fourier Transform Spectrometer (GIFTS) Sensor Module (SM) Engineering Demon- 
stration Unit (EDU) is a high resolution spectral imager designed to measure infrared (IR) radiances using a 
Fourier transform spectrometer (FTS). The GIFTS instrument employs three focal plane arrays (FPAs), which 
gather measurements across the long-wave (685-1130 cm -1 ) IR (LWIR), short/mid-wave (1650-2250 cm -1 ) IR 
(SMWIR), and visible spectral bands. The raw interferogram measurements are radiometrically and spectrally 
calibrated to produce radiance spectra, which are further processed to obtain atmospheric profiles via retrieval 
algorithms. In this paper, we emphasize on the algorithmic components of the calibration process, which are 
divided in multiple blocks. Fig. 1 illustrates the basic data processing blocks for the GIFTS SM EDLT Level IB 
algorithms. In the raw data processing block, raw interferograms are modified based on the detector nonlinearity 
characteristics. A complex finite impulse response (FIR) filtering and decimation procedure is then performed 
on the interferograms. The next block of procedures includes the steps for calibrating a single interferogram. 
Prior to the radiometric calibration, the filtered and decimated complex interferograms must be phase-corrected 
so that their imaginary spectra contain only noise. The phase correction is performed on all three types of 
measurements for a series of scans, i.e., the ambient blackbody (ABB) reference, hot blackbody (HBB) reference, 
and scene scans. To reduce the random noise embedded in the reference measurements, a spectral smoothing 
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method can be applied to the phase-corrected blackbody reference spectra. Next, the calibrated ABB, HBB, 
and scene radiance spectra are computed based on the previous results. An assessment of the noise equivalent 
spectral radiance (NESR) is derived from the calibrated ABB and HBB references. Two additional correction 
algorithms are implemented to compensate for the fore-optics offsets and FPA off-axis effects. In the multi-pixel 
algorithm block, we introduce an efficient method of estimating the noise performance of the entire FPA, as well 
as a random pixel selection strategy that can be employed for generating a representative pixel sample set. In 
the FPA calibration block, the calibration algorithms introduced previously in the second block are performed 
on the entire FPA. The interferogram data at various stages are retrieved, linked, and stored via the storage and 
management system. Detailed descriptions for each data processing block are given in the following sections. 
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Figure 1. The basic GIFTS SM EDU data processing blocks. 


2. RAW DATA PROCESSING 

In this section, we introduce the necessary procedures to be performed on the measured raw interferograms prior 
to the calibration. We will discuss how the detector nonlinearity affects the raw interferogram data as well as the 
method used to characterize and correct the distortion. We then continue onto the topic of complexing filtering 
and decimation techniques, which have been applied to the interferogram data. 

2.1 Detector Nonlinearity Correction 

Ideally, a given input radiance function I 0 (n ) and its corresponding measured output interferogram I(n) are 
linearly proportional. However, in practice, this relationship becomes less ideal due to the distortion caused by 
the detector nonlinearity effect. The input-output function can be modeled as a P th - order polynomial 

p 

y « F(x) = J2a p x p . (1) 

p=0 

The polynomial function in Eq. (1) and its coefficients can be obtained from a set of known data points on the 
input-output curve, i.e., (aq, j/i), (aq, 2/2), (xjv-i, yN-i), (a,’jv, Vn)- Fig. 2 depicts the effect of detector 
nonlinearity on the input I 0 (n). A sample polynomial curve with corresponding data pairs are also shown in 
the graph. The true amplitude of the input I 0 (n) can be estimated from the measured output I(n) by inverting 
the polynomial function or by interpolation. To characterize the GIFTS detector nonlinearity, a series of small 















Figure 2. An illustration of the input-output distortion caused by the detector nonlinearity effect. 



Figure 3. The detector nonlinearity characterization. 


signal linearity measurements were made at several blackbody temperatures. At each given temperature, the 
output response data is collected. After injecting a stim source AC signal into the blackbody signal, the output 
response is re- measured. This pair of measurements gives us the quantity of Ay = y on — y 0 ff- The variations 
in Ay indicate nonlinear behaviors in the detector response function. However, we do not have any information 
on the corresponding input values other than the assumption that the stim source is constant. Consequently, an 
extra procedure is necessary before we can make the nonlinearity correction. The input levels must be estimated 
in order to construct the input-output curve. Fig. 3 shows the stim on and off data pairs with the ideal linear 
and estimated nonlinear input-output curves. The algorithm for estimating the nonlinearity curve is summarized 
in the following steps: 


• Search for the linear region: 

- Compute output small and large signals for all pairs of measurements. 


- Compute the differences between adjacent elements from the previous results. 

- Compute the ratio r n . 

- The smallest value, r^, in r n indicates that the smallest change in Ay, which implies that the most 
linear portion of the curve is located between n = k and n = k + 1. 

- Estimate Ax. 

• Linear region estimation: 

- Estimate £ 0 //,n and x oritn for n = k, k + 1. 

• Forward region estimation: 

- For every pair of data above the linear region, we estimate their values from the previous pair using 
the finite difference approximation approach. 

Only one intercept point is constructed (one local extremum) to ensure the smoothness of the curve. 
The intercept point is placed at the midpoint between two adjacent pairs of data. 

• Backward region estimation: 

- Reverse the procedure for data pairs below the linear region. 


Fig. 4 depicts three regions of the estimated curve. Once the detector nonlinear response function is constructed, 
the true input interferogram data can be found using interpolations. 
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Figure 4. The estimation of the detector nonlinearity function. 


2.2 Decimation Filtering 

The next step involved in the raw data processing block is the procedure of complex filtering and decimation. 
For the GIFTS SM EDU data, we designed a linear-phase even-length complex finite impulse response (FIR) 
filter that has the form of 

h(n) = h r (n) + jhi(n), n = 0,...,N — l. (2) 

In the preceding expression, h r {n ) contains the real part of the filter coefficients, which is the impulse response 
of an optimum equiripple linear-phase FIR filter. 1 The imaginary part /q(n) is obtained by taking the Hilbert 
transform of h r [n). The output of the filter is complex in the interferogram domain, and does not contain 
an image spectrum in the frequency domain. The filter design parameters include: S p - passband ripple, S s - 


stopband attenuation, u> p - passband cutoff frequency, u> s - stopband cutoff frequency, and N- filter length. The 
Fourier transform of Ecp (2) can be written as 

= H r (e j “) + jH 2 (e^), (3) 

where 

H r (e j “) = F{h r {n)} and Hi(e> u ) = T{hi{n)}. (4) 

Next, we derive the imaginary impulse response hi(n) by taking the Hilbert transform of h r (n), for instance, 

hi(n) = TL{h r (n)} = h r {n) * hh(n). (5) 

The real and imaginary spectra can be expressed as 


H r {e> u ) = \[H{e iu ) + H*{e~ ju )} and 
j Hi(en = \[H(en~H*{e~n]. 

From the preceding expression, the imaginary spectrum also can be written as 

0 < w < 7T, 


or 


—7 r < uj < 0 , 

= H h {e j “)H r {e h “), 


where 

[j, -TV < UJ < 0. 

Note that Hh(e^ u ) is the spectrum of a Hilbert transformer, which can be seen as a 90-degree phase shifter. 


( 6 ) 

(7) 

( 8 ) 
(9) 

( 10 ) 


3. SINGLE-PIXEL CALIBRATION 

For a single pixel element on the FPA, the calibration procedures involve the following steps: phase correc- 
tion, spectral smoothing, responsivity calculation, radiometric calibration, NESR assessments, fore-optics offsets 
adjustment, and off-axis effects correction. 

3.1 Phase Correction 

The phase correction algorithm is applied to the decimated and filtered complex interferogram I m (n). Its complex 
spectrum N m (a) has the form of 

N m (a)=N(a)e j ^ +n(a), (11) 

where n(a) represents the noise vector, and N(c r) is the real and noise-free spectrum. To estimate the phase 
function </>(ct), we adopt the Forman-Vanesse-Steel (FVS) method. 2 In this technique, a Hamming window 
w(n), of length W, is applied to the complex interferogram with respect to the zero path difference (ZPD), i.e. , 
n = n zp d ■ The windowed interferogram is written as I w (n ) = w(n)I m (n), where 


w-i 


_ J °- 54 - 0.46 cos (j^rzy), \n-n zpd \< 2 , 
I 0, otherwise. 


( 12 ) 


The resulting interferogram I w (n) is circular shifted to obtain I' w (n) such that I' w { 0) = I w (n zp d)- The spectral 
response of I' w (n) is computed via the Fourier transform, i.e., = J-{I' w (n )}, from which, the phase 

function can be estimated as (j){a ) = tan -1 . By applying the phase estimate, the phase-corrected 

spectrum becomes 

N{a) = N(a)e jWa) -^ a)) + n{a)e~ j ^ a) , (13) 

and its resulting real and imaginary components are given as 

Re{A(cr)} « N(a) + Re{n(cr)} cos(0(cr)), (14) 

Im{A(cr)} « Im{?r(a)} sin(^(cr)). (15) 

Note that the imaginary part of the spectrum contains only noise, which can be discarded. 



3.2 Spectral Smoothing 

The phase correction is performed on three types of spectra, i.e. , ABB, HBB, and scene. An optional smoothing 
algorithm can be applied to the blackbody references to improve the noise performance of the calibrated radiances. 
Two smoothing techniques are being considered for this purpose: the least-squares smoothing filter and the 
moving average (MA) filter. The MA smoothing method approximates the function within the predefined 
window as a constant whereas the least-squares algorithm fits a polynomial to all points within the window. The 
MA method does not preserve higher order moments, however, it is less demanding computationally. See 1-3 for 
detailed reviews of these two techniques. When the final noise performance of these two methods are compared, no 
significant differences were observed; therefore, the MA smoothing technique is chosen due to its simplicity. Fig. 5 
(a) illustrates the phase corrected and smoothed longwave spectra for the real and imaginary blackbody references 
at ambient and hot temperatures, and Fig. 5 (b) shows the phase response functions after the correction. 
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LW Phase Response (Corrected) Pixel (64, 65), 25 scans, ABB00494_HBB00483 ,res.=0.57cm"' 
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Figure 5. Phase corrected longwave blackbody spectra: (a) real and imaginary blackbody spectra at ambient (260K) and 
hot (286K) temperatures, and (b) phase response functions. 


3.3 Radiometric Calibration 

To obtain the calibrated radiances, the first step involves the calculation of the spectral responsivity, which can 
be computed from 


R( N H (a)-N A (a) 
[<T) B H {a)-B A {a)' 


(16) 


N H (a) = ^ Ef =1 W/^ct) and N A (a) = ^ A^cr) represent the mean HBB and ABB spectra averaged over 
S scans, respectively. Using the relation of 


N H (a)=R(a)(B H (a) + 0(a)), 


(17) 


the offset term can be solved from 


O(cr) 


N A (a)B H (a) - N H {a)B A {a) 
N H (a) - N A {a) 


(18) 


Next, the calibrated HBB, ABB, and scene radiance spectra can be determined, based on the responsivity and 
offset results, from 


Hf’V) 


Nf(*) 


N z H ’ A (a) 

R(a) 

Nf(a) 

R(a) 


- 0{c r), 
O(a). 


(19) 


( 20 ) 


respectively. The variable i denotes the scan number, i.e. , i = 1, 2, . . . , S. 


3.4 Noise Equivalent Spectral Radiance (NESR) 


This section describes the estimation of the noise equivalent spectral radiance (NESR) from calibrated blackbody 
references. The NESR is generally considered as a measure of the instrument noise performance, and can be 
estimated as the standard deviation of calibrated blackbody radiance spectra from multiple scans. The NESR 
definition is given by 


NESR H ’ A (a) 



E (*f ’ V) 


i=l 


B H ’ A (a)) 2 . 


( 21 ) 


Fig. 6 depicts the smoothed HBB and ABB NESR estimates over their corresponding non-smootlied NESRs 
computed from 25 interferogram scans. 
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Figure 6. Smoothed and non-smoothed longwave HBB and ABB NESR estimates computed from 25 interferogram scans. 


3.5 Fore-optics Offsets 

The previous calibration results are computed using internal blackbody references. To compensate for the offsets 
generated from fore-optics, a correction scheme is incorporated into the previous model. This can be achieved by 

jj a 

using additional data measurements collected from an extended source. Given C e (a) and C e (a) as calibrated 
extended blackbody sources averaged over all scans, we have 

C*' A (a) = R e (a)B?’ A (a) + O e (a). (22) 


By solving for R e (cr) and O e (a), the actual scene radiances can be approximated as 

EV) - °e(cr) 


Nf(a) = 


Re(cr) 


(23) 


Fig. 7 shows the calibrated longwave HBB, ABB, and scene radiance spectra after performing the fore-optics 
offsets correction. 
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Figure 7. Calibrated longwave HBB, ABB, and scene radiance spectra after the fore-optics offsets correction (25 scans). 

3.6 Off-axis Effects 

In this section, the correction of the FPA off-axis effects is presented. Since the GIFTS FPA contains 128 x 128 
detector elements, the distortion caused by the off-axis effects is too great to be neglected. The off-axis pixels are 
sampled at slightly shorter optical path differences (OPDs), which cause the spectra of these pixels to expand 
to slightly higher wavenumbers. 4 The correction can be formulated as a fractional sampling rate conversion 
problem, which can be solved via sine interpolations, i.e. , 


M — 1 

N i (o') = Sine (cr - a'), a' = (24) 

m = 0 ^ 

where cr is the on-axis wavenumber scale, a' is the off-axis scale, and / is the off-axis factor that is computed 
from the FPA geometry. However, the sine interpolation calculation is computational intensive, furthermore, 
the truncation window associated with the sine kernels may cause overshoots in the resulting data. It has been 
shown that the zero-padding interpolation in the interferogram domain is equivalent to the sine interpolation 
in the spectral domain except that the zero-padding method cannot evaluate the interpolation output at an 
arbitrary point. 5 To overcome this limitation, we designed the method of “over-padding”, in which, we assign 
an over-padding factor g that is closely correlated to the final resolution of the corrected spectra. For instance, 
if / = 0.9977, and the length of the spectrum is N, then the zero-padded interferogram length is the roundoff 
value of g * N/f. If g = 100, the actual off-axis factor of 0.997702 is obtained; if g = 1, then the actual value 
of / becomes 0.99758 due to roundoff errors. The corrected spectrum can be obtained by downsampling the 
Fourier transform of the over-padded interferogram by the factor g. Fig. 8 (a) depicts the spectra for two pixels 
before the off-axis effects correction, one pixel is located near the center of the FPA whereas the other is situated 
near the corner of the FPA. A better wavenumber alignment was achieved after performing the off-axis effects 
correction as shown in Fig. 8 (b). 


4. MULTI-PIXEL ALGORITHMS 

We have discussed the algorithmic procedures required for a single pixel measurement. To evaluate the entire 
FPA without completing all of the calculation steps, we designed a pixel inventory method, which estimates the 
FPA responsivity and noise distributions. 
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(a) (b) 

Figure 8. An illustration of the off-axis effects: (a) calibrated spectra for Pixel (1,1) and Pixel (64, 65) before performing 
the off-axis effects correction, and (b) after the off-axis effects correction. 


4.1 Pixel Inventory 


The responsivity for any detector element can be approximated from the normalized intensity value at the ZPD, 
which is given by 


R x ’ y = 


Im V { n z P d) 




V x,y , 


(25) 


where x , y are pixel numbers. The FPA responsivity estimate is shown in Fig. 9 (a). The FPA responsivity 
distribution within a predefined range of 0.8 < R x,v < 1.2 is plotted in Fig. 9 (b). To assess the FPA noise, each 
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Figure 9. Longwave pixel inventory: (a) FPA responsivity estimates, and (b) FPA responsivity distribution plot. 


interferogram is first normalized by its value at n = n zp d , i.e., 


Next, the estimated noise is computed from the root-mean-square value of the last L samples in each interferogram 
that is represented by 


N x ’ v = 


N 


1 M —1 

7 ^2 \l™ v (n)\ 2 V x,y. 


Similarly, the FPA noise and its distribution are shown in Figs. 10 (a) and 10 (b), respectively. 


(27) 
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Estimated Noise (Normalized by ZPD) Distribution for lwp_ABB00494_scan_01 


mean = 0.00036859 
median = 0.00034939 
stdev. = 7.8007e-005 
number = 12237 
num < 0.0001 = 8 
num > 0.0007 = 4139 
total outliers = 41 47 
% in range = 74.6887 
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Figure 10. Longwave pixel inventory: (a) FPA noise estimates, and (b) FPA noise distribution plot. 


4.2 Random Pixel Set Generation 

Once we have the knowledge of the FPA responsivity and noise estimates, a set of selection criteria is established 
by defining an acceptable range for the responsivity and noise. Thus, a random sample set of pixels can be chosen 
such that it reflects the statistical properties of the entire FPA. We can then perform the actual calibration 
procedures on this sample set. Fig. 11 (a) shows a sample longwave pixel selection scheme, in which 64 random 
pixels are selected with 4 pixels per tap. The selection is uniformly distributed and unique within each tap. The 
NESR estimates of this sample pixel set can be seen in Fig. 11 (b). 

5. FOCAL PLANE ARRAY (FPA) CALIBRATION 

The calibration process for an entire FPA can be summarized into the following steps: 

• Perform detector nonlinearity corrections for all pixels. 

• Implement the complex filtering and decimation process for all pixels. 

• Generate calibrated HBB and ABB data cubes from the extended blackbody source measurements. 

• Calibrate the scene data cubes using the internal blackbody references and results obtained from the 
previous step. 

• Compensate for the fore-optics offsets and off-axis effects. 

Fig. 12 (a) shows a particular frame of calibrated radiances before applying the off-axis correction. The circular 
pattern illustrates the misalignment caused by the off-axis effects. Fig. 12 (b) demonstrates an improved spectral 
alignment after the correction. 




2 — 
1 .8 - 
1.6 - 
^ 1.4 - 

\ 12 ' 
\ 1 - 

E 0.8- 
§ 

E 

“ 0.6 - 
0.4 - 
0.2 - 


LW ABB NESRs: 64 pixels, 25 scans, ABB00494 


700 750 800 850 900 950 1000 1050 1100 1150 

Wavenumber [cm -1 ] 

(b) 


Figure 11. A sample longwave random pixel set: (a) blue pixels: out-of-range noise; green pixels: out-of-range responsivity; 
black pixels: out-of-range noise and responsivity; gray pixels: satisfy both noise and responsivity constraints; and white 
pixels: random pixel selection, and (b) NESR estimation results of the sample pixel set. 



Figure 12. A calibrated longwave FPA frame: (a) before the off-axis effects correction, and (b) after the off-axis correction. 


6. CONCLUSION 

This paper describes Level IB algorithms for calibrating the GIFTS SM EDU data spectrally and radiomet- 
rically. Four basic processing blocks were introduced. Initially, the raw interferograms must be corrected for 
the distortion caused by the detector nonlinearity, then the complex filtering and decimation procedure was 
performed. In the next block, the basic calibration steps were carried out. These include the procedures for 
phase correction, spectral smoothing, radiometric calibration, NESR computation, correction algorithms for the 
fore-optics offsets and the off-axis effects. The next two blocks of algorithms are concerned with multi-pixel 
processing techniques. Methods for estimating the FPA’s performance, generating random pixel sample sets, 
and calibrating the complete FPA were discussed. 
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